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No spin diffusion in the spin 1/2 XXZ chain at T = oo 
"S": Numerical asymptotics 
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Abstract 

We analyze the recent numerical computations made by Fabricius, Low and Stolze 
to show that the long time behavior of the zz correlation function of the spin 1/2 XXZ 
chain at T = oo is very well fit by the formula t~ d [A + Be -1 ^ 1-1 ^ cos 0(t — to)] where d is 
substantially greater than 1/2. This confirms the conclusion that there is no spin diffusion 
in this model. 

1. Introduction 

The spin 1/2 XXZ chain of N site with periodic boundary conditions specified by the 
Hamiltonian 

1 - 

H = 2 £ (Wi + + A «+i) ( L1 ) 

i=l 

where a\ is the y = x,y,z Pauli spin matrix at site % is well known to be an integrable 
system. As such it is to be expected that all the spin correlation functions should be 
analytically computable. In the past decade there has been much work done to fulfill this 
promise but much still remains to be done. In particular there has been no published work 
done on the time dependent autocorrelation function at infinite temperature 

S(t; A) =< a z (t)a z (0) >= lim N ^ 00 2- N Tre- UH a z e ltH a z (1.2) 

other than the old result M that at A = 

S(t-0) = [J (2t)} 2 (1.3) 

where Jo(2t) is the Bessel function of order zero. 

On the other hand this infinite temperature correlation function is of great interest 
because of its relation to the theory of spin diffusion @- [§J which says that in the limit 
t — > oo if spin diffusion is present in the system then the asymptotic behavior of S(t, A) 
should be 

S(t; A) ~ At~ 1/2 . (1.4) 

In the absence of exact results this asymptotic behavior has been studied by numerical 
methods for well over 20 years 0-0 • By far the most accurate of these numerical studies is 
the recent work of Fabricius, Low, and Stolze who have published a beautiful high pre- 
cision study of the time dependent two spin correlation functions at various temperatures 
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of the spin 1/2 XXZ antiferromagnet obtained by extrapolating to the thermodynamic 
limit a finite size exact diagonalization study of systems of size up to N=16. From this 
study the authors were able to give precise estimates of S(t; A) for A = 1, cos(.3tt) for the 
time interval < t < 4.95 and they concluded that no evidence of the asymptotic behavior 



could be seen. This supports the conjecture made by the present author [|10| and by 
others [[ll| that the integrability of the XXZ chain forbids the spin diffusion asymptotic 
behavior (|1.4j ) and agrees with the absence of diffusion at low temperature found by 
Narozhny Q for -1 < A < 1. 

However, in || there is no positive statement given as to what the true asymptotic 
behavior of their data might be. Here we show that the numerical results of || can be 
beautifully fit with the form 

fit- d, A, B, 7, O, t ) = t~ d [A + Be- l{t - to) cos tt(t - t )\ (1.5) 

which was first suggested in || with d = 1/2. 



To test the validity of the form (|1.5|) we have computed 



g(t;d,A) = t d S(t;A)-A, (1.6) 

where S(t; A) is given by the numerical values of which were generously sent to us in 
electronic form, and compared the result with 

h(t;B,<y,tt,t ) = Be- l{t ~ to) cos Q(t -t ). (1.7) 

The best results of these comparisons are given graphically in Figs.l and 2 for A = 1 and 
cos(.37r). From these comparisons we conclude that for 2.2 <t< 4.95 the equality 

S(t;A) = f(t;d,A,B, 1 ,n,t ) (1.8) 

holds to a very high degree of precision where 

1) for A = 1 

d=.698, A = .245, B = .0581, 7 = .70, = 4.36, t = 1.90 (1.9) 

2) for A = cos(.3tt) = .587- •• 

d=.838, A = .21, S = .114, 7 = .354, = 4.06, t = 1.96 (1.10) 
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These results are to be compared with the exact asymptotics of A 
(0) as 



S(t;0) 



1 9/ 7I\ 1 

cos 2 (2t - -) = 

nt y A' 2nt 



1 + cos(4i 



7T , 



obtained from 



;i.n) 



which is of the form (|1.5| ) with 
3) A = 



d=l, A = B = -!- = .159---, 7 = 0, = 4, t = = 1.963 •• • (1.12) 

Z7T 8 

which is plotted in Fig. 3. 

In an attempt to provide some estimate of error in the fitting parameters we plot in 
Fig. 4 a fit of the case A = cos(.37r) obtained with 

d = .784, A = .196, B = .104 7 = .342, ft = 4.06, t = 1.96. (1.13) 

From comparing the fits of Fig. 2 and 4 we conclude that the estimates of Q and to are 
very stable while the exponent d varied by .045. However because the fit which better fits 
the points with larger t gives the larger value of t we feel that the conclusion that for 
t — > 00 the exponent d for A = cos(.37r) is definitely greater the value it has for A = 1 is 
certainly justified. 



2. Discussion 

The exponent d in ( p..9| )-( |TTT2 ) is always substantially greater that 1/2. Thus, 
if /(£; d, A, B, 7, f2, to) of ( J1.5| ) does indeed represent the true asymptotic behavior of 
S(t; A) then the spin diffusion form ( |1.4| ) is definitely eliminated. This is in agreement 
with and is the basis for the claim in the title of this paper. However, it must 
always be kept in mind that no study for a finite time interval can definitively claim to 
have seen the true t — > 00 behavior of a function. This obvious general statement is 
of particular importance here because numerical studies on the (nonintegrable) classical 
(S = 00) Heisenberg magnet |13[]-[]16| show that while a form like ( |1.5| ) held for times up 
to the order of 50 with an exponent d definitely greater than 1/2 that when times up to 100 
were considered the exponent d eventually became the spin diffusion value d = 1/2. The 
belief embodied in the conjecture of the absence of spin diffusion for an integrable model is 



that there will never be a scale at which the spin diffusion form (|1.4| ) sets in. We further 
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conjecture for the (nonintegrable) XXZ chains of arbitrary spin S > 1 that for suitable 
small times the form ( |1 . 5| ) will be an excellent approximation to the autocorrelation 
function but that there will exist some time scale such that for times larger than this scale 
spin diffusion will set in. This scale should increase as the spin S decreases and become 
infinite for the integrable case of S = 1/2. A similar remark holds for the addition of further 
nearest neighbor interactions to the XXZ chain which also destroy the integrability. 

We conclude by remarking that since the data of for T = oo is so very well fit by 
( |1.5| ) a similar simple fitting function must exist for finite T. In addition further numerical 
study should be able to determine the parameters in ( |1 . 5| ) as functions of A in the entire 
range < A < oo. As a guide to such a study we suggest that the parameter d(A) of 
( p. . 5|) is a monotonic nonincreasing function of A and that limA^oo^(A) = 1/2. 
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Fig. 1 Comparison of the numerical data of || for A = 1 with the form ( |1.5|) with 
d = .698, A = .245, B = .0581, 7 = .7, O = 4.36 ••• , t = 1.90. The points are 
the function g(t;,d,A) ( fOl ) obtained from the data of Pj. The continuous curve is 
h(t; B, 7, O, t ) of O with B = .0581, 7 = .70, O = 4.36, £ = 1.90. 
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Fig. 2 Comparison of the numerical data of |§ for A = cos(.37r) with the form 
(|L5|) with d = .838, A = .210, B = .114, 7 = .354, O = 4.06, t = 1.96. The points 
are the function g(t; , d, A) (|1.6| ) obtained from the data of ||. The continuous curve is 
h(t; B, 7, O, t ) of wi th B = .114, 7 = .354, O = 4.06, t = 1.96. 
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Fig. 3 Comparison of the exact result (|1.3|) for A = with the form ( |1 . 5|) with 
d = 1.0, A = B = i = .159 • • • , 7 = 0, O = 4.0, £ = ^ = 1-963 • • • . The points are the 
function g(t; , d, A) ( |1.6|) obtained from (pT3|). The continuous curve is h(t; B, 7, O, £q) of 
(PD with S = ^ = .159 • ■ ■ , 7 = 0.0, O = 4.0, t = ^ = 1.963 
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Fig. 4 Comparison of the numerical data of |§ for A = cos(.37r) with the form 
(|L5|) with d = .784, A = .196, B = .104, 7 = .342, O = 4.06, t = 1.96. The points 
are the function g(t; , d, A) (|1.6| ) obtained from the data of ||. The continuous curve is 
h(t; B, 7, O, t ) of O with B = .104, 7 = .342, O = 4.06, t = 1.96. 
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